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Abstract 

This paper compares temperature field 
measurements from selected experiments on a 
single row, and opposed rows, of jets 
injected into a ducted crossflow with 
profiles calculated using an empirical 
model based on assumed vertical profile 
similarity and superposition, and 
distributions calculated with a 3-D 
elliptic code using a standard K-E 
turbulence model. The empirical model 
predictions are very good within the range 
of the generating experiments, and the 
numerical model results, although 
exhibiting too little mixing, correctly 
describe the effects of the principal flow 
and geometric variables. 


Introduction 

The problem of jets-in-crossf low has been 
rather extensively treated in the 
literature, to the point that it can almost 
be called a "classical" three-dimensional 
flow problem. Although these studies, of 
chimney plumes, exhaust from V/STOL 
aircraft, discharges into rivers and 
streams, film cooling, and hot gas dilution 
etc., have all contributed additional 
understanding of the general problem, the 
information obtained in any given study is 
determined by its motivating application, 
and may not satisfy the specific needs of 
other diverse problems. 

One factor making the combustor dilution 
zone jet-in-crossflow application unique is 
that it is a confined mixing problem, with 
from 10 to 50 percent of the total flow 
entering through the dilution jets. The 
result is that the equi librium temperature 
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of the exiting flow may differ 
significantly from that of the entering 
mainstream flow. To control or tailor the 
combustor exit temperature pattern it is 
necessary to be able to characterize the 
exit distribution in terms of the upstream 
flow and geometric variables. This requires 
that the entire flow field be either known 
or modeled. , . 

Considerations of dilution zone mixing in . 
gas turbine combustion chambers have 
motivated several studies of the mixing 
characteristics of single and multiple jets 
injected into a crossflow in a constant 
area duct.r*«* 1-t? Recently,, 
experimental, results have been reported 
which provide insight into the. effect' of 
several geometric and flow variations . ... 

characteristic of most gas turbine ,, ; . 
combustion chambers, namely a variable. ... 
temperature mainstream, flow area 
convergence, and opposed rows of jets, 
either in-line or staggered. 10-1 ** 

From the data of references 1, 11, and 14, 
an empirical model was 
developed'*'® and 

extended 11 ' 1 '* for predicting the 
temperature field downstream of a single 
row, or opposed rows, of jets mixing with a 
confined crossflow. An interactive 
microcomputer program (Apple D0S : 3-3) based 
on the model of reference 5 was used in 
reference 9 to study the effects of . 
separately varying the momentum flux ratio, 
density ratio, orifice size, and orifice 
spacing, and to identify the relationship 
among these parameters which 
optimizes" the mixing. 

Although empirical correlation of 
experimental data can provide an excellent 
predictive capability within the parameter 
range of the generating experiments, they 
must be used with caution, or not . at all, 
outside this range. Physical modeling, in 
various levels of soph i sti cation 7 and 
complexity, may be used to obviate. this 
weakness. In this regard, several one and 


+ Although it is recognized that a uniform 
temperature distribution may not always be 
desired, "optimum" is used here, as in 
refs. 2, 5, 9, & 13, to identify geometries 
which, for a given flow condition, result 
in a uniform temperature distribution in a 
minimum downstream distance. 
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two dimensional integral and differential 
jet— in— crossflow models are extant which 
have been shown to give, for example, 
trajectory predict! ons that are in good 
agreement with experiments. These models 
may provide insight into the dominant 
physical mechanism <s) , and predict some of 
the characteristic parameters well, but 
they rarely provide sufficient information 
to quantify all of the important flow field 
variables (e.g. temperature, pressure, 
velocity, and turbulence quantities) in 
three coordinate directions. 

Recently, rapid advances have been made in 
the capability of computational fluid 
dynamics models and their application to 
complex flows such as 
jet <s> -in-crossflow. ie5 - 1Q These 
models are, however, still in the 
development and verification stage. They 
have been shown to be capable of predicting 
trends in complex flows, but their 
capability to provide accurate, 
quantitative and grid independent 
calculations of these flows has not yet 
been demonstrated. 

The present paper will compare temperature 
field measurements from selected cases in 
references 11 and 14 with distributions 
calculated using an empirical model based 
on assumed vertical profile similarity and 
superposition, 1 1 - and a 3— D 
elliptic code with a standard K-E 
turbulence model, 17 The results 
will show the capability (or lack thereof) 
of the models to predict the effects of the 
principal flow and geometric variables. 


flow area convergence. 1 1 * lJ * A 
major weakness of these (and previous) 
correlations is that their form precludes 
their use for semi confined flows (large 
Ho/D or S/D), single jet flows, or 
flows in which it is known a priori that 
the primary assumptions in the model will 
be violated. 

Numerical models do not have these 
limitations, and in addition provide 
calculations for all flow quantities of 
interest, not just those which happen to 
have been empirically correlated. The code 
used in this investigation is based on the 
USARTL 3— D model, and uses 
pressure and velocities as the main 
hydrodynamic variables. This code, or 
similar versions thereof, has been used in 
previous validation and assessment 
studies. 1S3 ~ 10 

The governing equations are represented by 
finite difference approximations on a 
staggered grid system. The differencing 
technique employed is hybrid for convective 
terms with central differencing of all 
other terms. The vel oc it y— pressure coupling 
is handled by the SIMPLE algorithm of 
Patankar and Spalding. 20 * 31 
Uniform velocities, and mass flow rates 
were used at all in— flow boundaries. The 
turbulent Schmidt Number was 0.9, the RMS 
turbulence intensity was chosen to be 5.5 
percent of the local mean velocity, and the 
inlet length scale was 2 percent of the jet 
diameter and duct height for the jet and . 
mainstream respectively. Standard values of 
the turbulence constants were used in all 
calculations. 


Description of Models 

The empirical model for the temperature 
field downstream of jets mixing with a 
confined crossflow is based on the 
observation that properly 

non— dimensional i zed vertical temperature 
profiles everywhere in the flow field can 
be expressed in the following self— similar 
form® 


(T mm x ^ ~~ T ) 

(T m .H^- - T<=) 

- <ln2) (y-v^ ) 3 
= exp (« 1 / 2- >/ -)2 


(1) 


where T is the temperature at vertical 

location y, and 

Wi/ 2 >/- , T c , and 

y c are scaling parameters (see 

Nomenclature). Correlations have been 

developed for each of these in terms of the 

independent variables J, S/D, 

Ho/D, Z/S, and X/Ho. The 

correlation of reference 5 for a single row 
of jets in a uniform temperature crossflow 
have been extended for predicting the 
temperature field downstream of single or 
opposed rows of jets, either in-line or 
staggered, injected into an isothermal or 
non— i sothermal mainstream, with and without 


Numerous uncertainties are currently 
present in 3— D calculations due to, for 
example, numerical diffusion, unmeasured 
(and hence assumed) boundary conditions, 
and turbulence model assumptions. The 
results shown here are not intended to 
represent the "best" agreement possible 
from numerical models at this time, as 
better temperature field agreement could 
undoubtedly have been achieved by adjusting 
model constants and/or inlet boundary 
conditions. Since this was not necessary to 
satisfy the present objective of evaluating 
the potential of these codes vis— a— vis 
combustor dilution zone flowfields, and due 
to the substantial uncertainties in the 
numerical calculation and because the mean 
temperature was the only parameter 
compared, no adjustments were made. 

Results and Discussion 

Figure 1 shows a schematic of the dilution 
jet flowfield of interest in this paper. 

The temperature field results are presented 
in three-dimensional oblique views of the 
temperature difference ratio, THETA, where 

THETA = (Ja - T ) 

(T m - Tj ) 

( 2 ) 

Note that THETA = 1 if the local 
temperature is equal to the jet 
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temperature, and THETA = O if the local 
temperature is equal to the mainstream 
temperature. The equilibrium THETA for any 
configuration is equal to the fraction of 
the total flow entering through the 
dilution jets, Wj/w-r- 

In the 3-D plots the temperature 
distribution is shown in planes normal to 
the main flow direction. The coordinates y 
and z are, respectively, normal to and 
along the orifice row in this constant x 
plane. The orifice configurations 
investigated are shown in figure 2. For 
clarity and consistency of the visual 
presentation, the THETA distributions are 
shown over a 2S span in the z-di recti bn, 
with the plane between jets, the midplane, 
at the edge of the oblique plots. 

The following paragraphs, and the 3-D plots 
in figures 3 to 11, compare the 
experimental data with the empirical and 
numerical model results, in the context of 
the effects of the primary independent 
variables. The flow and geometry conditions 
corresponding to these figures are given in 
table 1. Isotherm contour plots and 
centerplane temperature profile plots for 
these conditions are given in reference 22. 

Variations with Orifice Size and 
Spacing. At constant orifice area, 
changes in orifice size and spacing can 
have a significant influence on the THETA 
profiles. This is shown in parts a) and b) 
of both figures 3 and 4, where jets from 
closely spaced small orifices 
under -penetrate and remain near the 
injection wall, and jets from widely spaced 
larger orifices over — penetrate and impinge 
on the opposi te wal 1 . 

The empirical model reproduces the data 
very well in the small orifice case, since 
the data are consistent with the major 
assumption in the empirical model, that all 
vertical temperature distributions can be 
reduced to similar Gaussian profiles. The 
empirical model does not do as well in the 
larger orifice case however, as the 
impingment of the jets on the opposite wall 
results in vertical profiles which are not 
similar. 

The numerical model calculations made with 
approximately 20,000 nodes, although in 
qualitative agreement with the data, show 
temperature gradients that are too steep, 
especially in the transverse direction. 
Under — prediction of the mixing was seen in 
the single— jet calculations of reference 15 
also, where it is shown that the isotropic 
turbulence model under —est i mates the 
intensity. The result in figures 3 and 4 is 
typical of the numerical model calculations 
to be shown in this paper. 

For the smal 1— ori f ice case a coarse— grid 
calculation using less than 6000 nodes was 
also performed. The numerical results in 
figures 3a) and 4a) illustrate the 
significant influence of grid selection on 
the solution obtained, and the smearing of 


the profiles which can occur as a result of . 
numerical diffusion. 


Coupled Spacing and Momentum Flux 
Ratio. Examination of the experimental 
data revealed that similar jet penetration 
is obtained over a range of momentum flux 
ratios, independent of orifice diameter, if 
orifice spacing and momentum flux ratio are 
correctly coupled. 2 * 0 * 9 * 13 For 
example, low momentum flux ratios require 
large, widely spaced holes, whereas smaller 
closely spaced holes are appropriate for 
high momentum flux ratios, as shown in 
figures 5 and 6. 


In general, jet penetration and centerplane 
profiles are similar when the spacing is 
inversely proportional to the square root 
of the momentum flux ratio, i.e. : 


S/Ho = C/< 60 RT ( J) ) 


<3) 


For single— side injection, the centerplane 
profiles are approximately centered across 
the duct height and. approach an isothermal 
distribution in the minimum downstream 
distance when C=2.5. This appears to be 
independent of orifice diameter, as shown' 
in both the calculated and experimental 
profiles in figures 7 and 8. 

Values of C in equation (3) which are a . 
factor of >2 smaller or larger than 
the optimum correspond to under-penetration 
or over — penetration respectively (e.g. 
figures 3 and 4). 


In all of the combinations shown in figures 
5, 6, 7, and 8 the empirical model results 
are in very good agreement with the data, 
as the data. are consistent with the 
Gaussian profile assumption. The numerical 
model calculations using approximately 
20,000 nodes for these cases show jet 
penetration which is in good agreement with 
the data,, but the mixing is otherwise* »« • 
under — predicted as described previously.. 

Examination of the empirical, model results : 
in figures 3 to 8 in the context of 
equation (3) suggests that in .general the 
empirical model would be expected to 
provide good temperature field predictions 
for single-side injection when 
1< (S/Ho) (sonTU) K5. It 

is significant to note that the numerical 
model is not subject to this inherent 
limitation. 


Variable Temperature Mainstream. The 
influence of a non-i sother mal mainstream 
flow on the profiles for medium momentum 
flux ratios with S/H 0 =-5 and 
H 0 /D=4 can be seen by comparing 
figure 9 with figures 7b) and 8b) - The 
shape of the experimental profiles in 
figure 9 suggests modeling them as a 
superposition of the upstream profile and 
the corresponding jets-in— an— isothermal 
mainstream distribution. 13 The 
hottest temperature in the mainstream flow 
was used as T m in the definition 
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of THETA (see equation (2)) for this case. 

This gives a good approximation, as seen in 
the empirical model results. The agreement 
is only first-order however, since with a 
variable temperature mainstream there can 
be cross-stream thermal transport due to 
the flow of mainstream fluid around the 
jets ( and hence to different y locations), 
and this is not accounted for in 
superimposing the distributions. This 
becomes apparent if the local mainstream 
temperature, T m (x,y), is used in 
the definition of THETA in equation (2). 

In the variable temperature mainstream case 
the numerical model results agree well with 
the experimental data, especially on the 
jet centerplane, but the transverse mixing 
is underpredicted, as in the corresponding 
isothermal mainstream case in part b) of 
figures 7 and 8. 

Opposed Rows of In-line Jets. For 
opposed rows of jets, with the orifice 
centerlines in-line, the optimum ratio of 
orifice spacing to duct height is one-half 
of the optimum value for single— side 
injection at the same momentum flux 
ratio- 13 As as example consider 
the single-side case with S/Ho-. 5 
and H 0 /D=4 in figures 7b) and 8b), 
and the opposed row of in-line jets with 
S/Ho—. 25 and H 0 /D=8 in 

figure 10. Note that the product of the 
spacing and the square root of the momentum 
flux ratio <C in equation (3)) is 1.25 for 
this case, or one-half of the optimum value 
for single— side injection. 

The empirical model predicts the 
opposed— jet case very well, verifying the 
primary assumption that the effect of a 
plane of symmetry is similar to that of an 
opposite wall. 3 -* 9 * 13 Note that the 
experimental profiles on both sides of the 
plane of symmetry support the Gaussian 
profile assumption. The numerical model 
results show the steep transverse and 
lateral gradients indicative of too little 
mixing, as seen in almost all of the 
previous calculations also, but the 
penetration and profile shape are in good 
agreement with the data. 

Opposed Rows of Staggered Jets. For 
opposed rows of jets, with the orifice 
centerlines staggered, the optimum ratio of 
orifice spacing to duct height is double 
the optimum value for single-side injection 
at the same momentum flux ratio. 13 
As an example consider the single— side case 
with S/H 0 =.5 in figures 7b) and 
8b) (C — (S/H 0 ) ( ganT ( J ) ) 

— 2. 16) , and the opposed row of staggered 
jets with S/H 0 =l in figure 11 <C = 

5.25) . 

The empirical model does not handle this 
complex case well, as the fluid dynamic 
interactions here are not amenable to a 
direct extension of the simple Gaussian 
profile and superposition type modeling 
appropriate for most of the single— side and 


opposed- jet cases of interest. The 
numerical model calculations, are not in 
appreciably better agreement with the data 
than the empirical model results, as the 
mixing is under — predicted here as in the 
previous cases. 

As in most of the previous numerical 
calculations, approximately 20,000 grid 
points were used, but note that twice as 
many grid nodes were required in the 
transverse direction in this case, and that 
the number of axial grid points was 
correspondi ngly reduced. The number of 
axial nodes used in this calculation is 
approx imately the same as that used in the 
coarse grid calculations shown in figures 
3a), 4a), 5b), and 6b). It follows that the 
numerical diffusion in the x-direction 
would be comparable between these 
calculations. 

The numerical result for the staggered jet 
case is encouraging in that this method is 
not restricted by profile assumptions as is 
the empirical model, and the numerical 
model results should improve with overall 
improvements in the capability of the 3— D 
codes. 


Summary of Results 

The present paper compares temperature 
field measurements from selected 
experiments, wherein a single row, or 
opposed rows, of jets were injected into a 
ducted crossflow, with distributions 
calculated with an empirical model based on 
assumed vertical profile similarity and 
superposition and with a 3— D elliptic code 
using a standard K— E turbulence model. The 
empirical model calculations of the 
temperature field are very good within the 
parameter range of the generating 
experiments whenever the primary 
assumptions in the model are supported by 
the experimental data. 

Although the numerical model calculations 
consistently exhibit too little mixing, the 
trends which result from variation of the 
independent flow and geometric variables 
are approximated correctly. Codes with 
improved numerics, accuracy, and turbulence 
models should provide more quantitive 
predictions. Of note is the fact that the 
3-D codes can provide calculations for 
complex flows that are outside the range of 
available experiments, or for which the 
assumptions in the empirical model are 
known a priori, to be invalid. Also 
important is the fact that the 3— D codes 
provide a prediction for all flowfield 
quantities (i.e. velocity, turbulence, 
etc.), not just those that have been 
empirically modeled. 
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Nomencl aturg 


t 


Cc* 

D 

DR 


Ho 


J 

H 

R 

S 

T 


T c 

T, 

Tm 


T 


THETA 

U 

U m 

V„ 

Wj/Wt 


X 


y 


z 


= jet-to— mainstream area ratio 

= (Pi /4) / ( <S/Ho) (Ho/D) 2 > Tor one— side injection 
= (Pi /2> / ( (S/H 0 > (Ho/D) 2 > Tor two— side injection 

— oriTice discharge coeTTicient 
= oriTice diameter 

= (D) (aawr (Cd ) ) 

= jet— to-main stream density ratio = (T m /Tj) 

= duct height 

= jet— to— mainstream momentum Tlux ratio = (DR)(R) 2 
= jet— to— mainstream mass Tlux ratio = (DR) (R) 

= jet -to— mainstream velocity ratio = (Vj/U m > 

— spacing between oriTice centers 
= temperature 

= temperature at y c 
= jet exit temperature 
= mainstream temperature 

= maximum temperature above ( + ) or below (— ) the centerline 

— T m — (Tm — Tj) (THETA m m' , ' / ); see reT . 5 
= (Tm - T)/(Tm - Tj) 

— velocity 

* mainstream velocity 
= jet velocity 

= jet— to— total mass Tlow ratio 

— W(<DR) <J))) <Cd) (A</A„) 

1 ♦ (sort (DR) (J) > ) (Cd ) (Aj /A m ) 

= jet halT-widths above (+) or below (-) the centerline; see reT. 5 

= downstream coordinate 
= O at injection plane 

= cross— stream (radial) coordinate 
= O at wall 

= ycr at location oT minimum temperature in a line x=const, z=const 

= lateral (circumTerential ) coordinate 
= O at centerplane 


i 


$ 
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Table 1. Flow and Se owe try Conditions 


Figure 

S/Ho 

Ho/D 

Aj /Am 

Co 

DR 

J 

Wj/%iT 

(S/Ho) (j 

3a, 

4a 

.25 

8. 

.05 

-60 

2. 1 

22.4 

.17 

1.18 

3b, 

4b 

1.0 

4. 

.05 

.67 

2.2 

23.5 

. 19 

4-85 

5a, 

6a 

1.0 

4. 

.05 

.73 

2.1 

5.3 

.11 

2.30 

5b, 

6b 

.25 

8. 

.05 

.61 

2.3 

92.7 

.30 

2.41 

7a, 

7b 

.5 

5.7 

.05 

-71 

2-2 

25.4 

.21 

2.52 

7b, 

8b 

.5 

4. 

- IO 

.61 

2. 1 

18.6 

.27 

2.16 

9 


.5 

4. 

- IO 

.61 

1.8 

31.3 

.31 

2.80 

10 


.25 

8. 

. IO 

.65 

2.1 

25.0 

.32 

1.25 

n 


1.0 

4. 

. IO 

-65 

2. 1 

27.6 

.33 

5.25 
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S/Hq: .5 1 .25 .5 

Hq/D: 4 4 8 5.7 

Figure 2. - Orifice configurations. 







Figure 4. - Effect of varying orifice spacing at constant area; Aj/Am ■ 0. 05, X/Hg 




NUMERICAL 

EXPERIMENTAL EMPIRICAL 40x23x21 







Figure 6. - Temperature distributions for coupled spacing and momentum flux ratio. Aj/Am * 0.05, X/Hg 





NUMERICAL 

EXPERIMENTAL EMPIRICAL 35x33x17 



Figure 7. - Effect of varying orifice diameter at constant spacing-, S/Hg - 0. 5, X/Hg = 0. 5. 






NUMERICAL 



(b) H 0 /D - 4, Aj/Am - 0. 10, J - 18. 6. 

Figure 8. - Effect of varying orifice diameter at constant spacing; S/H 0 - 0.5, X/H 0 





Figure 10. - Temperature distributions for opposed rows of in-line jets: S/Hq = 0. 25, Hq/D » 8, Aj/Am - 0. 10, J = 25. 0. 








NUMERICAL 

EXPERIMENTAL EMPIRICAL 22x27x33 
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